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' Abstract. We report the final analysis of a search for microlensing events in the direction of the Andromeda galaxy, 
which aimed to probe the MACHO composition of the M31 halo using data collected during the 1998-99 observational 
campaign at the MDM observatory. In a previous paper, we discussed the results from a first set of observations. Here, 
^ I we deal with the complete data set, and we take advantage of some INT observations in the 1999-2000 seasons. This 
merging of data sets taken by different instruments turns out to be very useful, the study of the longer baseline available 
, allowing us to test the uniqueness characteristic of microlensing events. As a result, all the candidate microlensing 

■ events previously reported turn out to be variable stars. We further discuss a selection based on different criteria, aimed 
', at the detection of short-duration events. We find three candidates whose positions are consistent with self-lensing 

■ events, although the available data do not allow us to conclude unambiguously that they are due to microlensing. 
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1. Introduction 

Gravitational microlensing is a powerful tool for the de- 
tection of dark matter in galactic haloes in form of MA- 
CHOs ( |Paczynski 1986| l . Intensive searches in our Galaxy 
have shown that up to 20% of the halo could be formed 
by objects of aroun d M ~ QAMq l|Alcock et al. 20001 
ILasserre et al. 20n?))l . This result is still debated (e.g. 
IJetzer et al. 2002)l . and remains to be confirmed. For this 
purpose it is useful to probe the MACHO distribution 
along different lines of sight. 

A survey of M31, nearby and similar to our own 
Galaxy, can supply several insights into this problem. 
Briefly, it tests a different line of sight in our Galaxy, 
it allows to probe M31's own halo globally and finally, 
the inclination of the M31 disk should give an unmis- 
takable signature of microlensing events IjCrotts 19921 
IBaillon et al. 19931 IJetzer 1994|l . 

Several collaborations have undertaken a search 
for microlensing events towards M31 in the past years: 
AGAPE HAnsari et al. 1999|l; MEGA IjCrotts et al. 20001 : 
POINT-AGAPE IIAuriere et al. 2001 |l : WeCapp 
IjRiffeser et al. 2001|l . In this same framework we 
have already reported on the analysis of a partial set of 
the MDM data in Calchi Novati et al. (2002, hereafter 
[Paper \\ . Furthermore, the POINT- AGAPE collaboration 
has given evidence of 4 short high S /N ratio microlensing 
candidates (jPauhn-Henriksson et al. 20021 2003). 

In the present paper we conclude the analysis of the 
MDM data started in [Paper I| We complement the MDM 
data with INT observations taken by the POINT- AGAPE 
collaboration. Furthermore, we discuss the issue of selec- 
tion criteria for microlensing events, especially with re- 
spect to the background of variable stars, and we present 
results from an analysis based on different criteria. In ^ 
we summarize the observational setup and give some de- 
tails of the data analysis, dealing in particular with some 
improvements in the selection procedure; in |3| we intro- 
duce the INT data and discuss how we use the longer 
baseline they give us to construct a strong supplementary 
selection criterion; f^is devoted to the results of the anal- 
ysis; in fjSlwe present our conclusions. 

2. Observations and data analysis 

2.1. Observational setup 

We analyse data collected on the 1.3 m McGraw-Hill Tele- 
scope, at the MDM observatory, Kitt Peak (USA), towards 
the two sides of M31 including the bulge^. Two fields are 
observed, located in a = OOh 43m 24s, 5 = 41°12'10" 
(J2000) ("Target", whose data have been the object of 
the analysis presented in [PaperU ) and a = OOh 42m 14s, 
5 = 4r24'20" (J2000) ("Control"). The Target field is 
centred on the far side of M31 while the Control field is 

^ The data are shared with the MEGA collaboration. 



centred on the near side. Throughout the observations a 
2048 X 2048 pixel CCD camera, with a field of view of 
17' X 17', has been used. 

Two filters, similar to standard R and / Cousins, have 
been used in order to test achromaticity. Furthermore, this 
colour information gives us the chance to have a better 
check on red variable stars, which can contaminate the 
search for microlensing events. 

We analyse data taken starting from October 1998 to 
the end of December 1999. While the baseline of the Con- 
trol field is the same as that of the Target field, the Control 
field has only about 20 nights of observations for both fil- 
ters^ , which is only about half as much data as the Target 
field. The average value of the seeing is ~ l."6. 

The photometric calibration is done with respect to 
a sample of reference secondaries identified in our frame 
( [Magnier et al. 19931 ). As the photometric conditions for 
the reference images in the two fields are similar, the zero 
point calibrations reported in [PaperTI still hold. 

2.2. Reduction and analysis 

The reduction procedure and the candidate selection were 
discussed in detail in |Paper"H Here we just recall the basic 
points and discuss in more detail some aspects that have 
since been improved. 

2.2.1. Pixel Lcnsing 

Due to the distance of the target, the potential sources 
of microlensing events are not resolved stars. We use 
the pixel lensing technique developed by the AGAPE 
group to detect flux variations of unresolved sources 
l|Ansari et al. 1997jl ^. In order to cope with photometric 
and seeing variations, we first choose for each filter a ref- 
erence image (the R reference image is also the geometric 
reference for both filters) and then calibrate the flux* of 
all other images with respect to the reference image by 
means of a linear correction. Our seeing correction is em- 
pirical and does not require the evaluation of the PSF of 
the image | |Paper l| ). 

2.2.2. Bump detection 

Following [Paper I| we look for bumps (in R band) by a 
statistical analysis of the light curve. To this end we con- 

Each night ~ 20(11) images are taken in the R{I) filter, 
which are then averaged to get the single image per night that 
we use in the following analysis. 

^ An alternative method based on image subtrac- 
tion is currently used by the MEGA collaboration 
( jTomaney fc Crotts 1996t . 

* With the notable exception of geometrical alignment, 
throughout the analysis we substitute for the pixel value the 
sum of the fluxes taken in a square of 5 x 5 pixels around the 
central pixel, i.e., the corresponding superpixel value. 
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struct the two estimators for the significance of a variation 
L and Q. We define 

L = -ln(njet„„pF($|$ > given ^bkg, , (2.1) 

where 



P(4>|$ > $j 
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$j and (jj are the flux and its associated error in a su- 
perpixel at time tj, ^tkg is an estimator of the basehne 
level, and a "bump" is defined as a variation with at least 
3 consecutive points more than 3 a above the baseline. We 
define 
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where Xconst is the calculated with respect to the 
constant-flux hypothesis and Xpacz is the calculated 
with respect to a Paczyhski fit. We stress that Q is evalu- 
ated along the entire light curve, while L is evaluated only 
along the bump. 

In _Paper 1| we carried out the analysis to detect clus- 
ters of variable light curves using L. We eventually re- 
tained from the sample of selected light curves (the ones 
with the highest value of L in a given cluster) only those 
characterized by Q > Qthresh- 

Here we follow a somewhat different approach in that 
we select clusters related to flux variations by demanding 
Q > Qthresh- Compared to the selection based on L, we 
are biased in favor of monobump-shape variations. More- 
over, we exclude all those light curves showing spurious 
variations for which a high value of L can happen to be 
induced, e.g., by an underestimation of the baseline level. 

To conclude with the bump detection, we evaluate for 
all the pixels in each cluster L, and we then proceed as 
in |Paper"T| This is a better approach because L is based 
on an evaluation of the baseline more appropriate to this 
purpose, that is, to selecting the light curve showing the 
maximum flux deviation from the baseline. 



2.2.3. Shape analysis 

For microlensing events the flux variation must be unique, 
it must follow (in the point-like and uniform-motion ap- 
proximations) the symmetric Paczyhski shape, and it 
must be achromatic. 

The flux as a function of time for an event with am- 
plification A(t) and an unlensed source (f>* can be written 
as 

^{t) = (bbka + [A{t)~l]^*, (2.4) 

where 4'bkg is the background fiux including (j)* . 

The achromatic shape analysis is carried out as in 
[Paper 1| We perform a 7-parameter Paczyhski fit to the 
two bands simultaneously (background level (pbkg, source 
flux without amplification 0*, and the parameters of the 



amplification: time of maximum amplification to, Einstein 
time ^E, and the impact parameter Umin)- In the absence 
of direct knowledge of the source flux from e.g., a Hubble 
Space Telescope (HST) image fe.g.. lAuriere et al. 2001|l or 
an indirect measurement of it from a high signal-to-noise 
ratio hghtcurve fe.g. lPaulin-Henriksson et al. 2003|l . there 
is a degeneracy among the three parameters (jf , t^, and 
Umin HGould 1996|l . For these cases we retain two "eff'ec- 
tive parameters", the time width = ^1/2 (^E, Umin) 
and the excess at maximum with respect to the back- 
ground ls.4>max = ^4>niax {4>* , Umin) (from which wc eval- 
uate the "magnitude at maximum" Rmax and the colour 
at maximum {R — I)c)- For an event to be achromatic 
we have to test whether, along the bump, the ratio of 
the deviation of the flux from the background in the two 
bands remains constant in time (A07?/A0/ — (j)*]i/4>*j — 
constant)'^. 

Furthermore, we carry out the Durbin- Watson test 
l|Durbin & Watson 1951|l . which is sensitive to time cor- 
relation among consecutive residuals. In effect, it is able 
to check whether consecutive points lie, for instance, all 
above or below the best Paczyhski fit shape, such be- 
haviour being characteristic of variable stars (and for 
which the test can say little). The test is based upon 
the evaluation of a Durbin- Watson coefficient, dw, which 
must lie in an interval that depends on the number of 
points along the light curve, e.g. for 20 points it must be 
1.41(1.15) < dw < 2.59(2.85) at the 10% (2%) level of 
significance. 



3. The baseline of INT data as a test for bump 
uniqueness 

3.1. The problem oj variable sources 

A main problem in the search for microlensing events is 
the estimate of the background noise given by variable 
sources. This is particularly serious in the case of pixel 
lensing for two main reasons. First, the class of stars to 
which we are in principle most sensitive are the red giants, 
for which a large fraction are variable stars (regular or 
irregular). Second, as we look for pixel flux variations, it 
is always possible to collect (in the same pixel) light from 
more than one source whose flux is varying. 

Thus, in the analysis we are faced with two problems: 
large-amplitude variable sources whose signal can mimic 
a microlensing signal, and variable sources of smaller am- 
plitude whose signal can give rise to non-Gaussian fluctu- 
ations superimposed on the background or on other phys- 
ical variations. 



^ Note that this same test can be applied to resolved-source 
microlensing, provided that "background" is replaced by "base- 
line", i.e. the combined flux from the true source plus any un- 
resolved blended companions. That is, blending breaks achro- 
maticity only for the ratio of the total fluxes in two bands, not 
differences in these fluxes from baseline. 
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In [Paper I| we followed a sufficiently conservative ap- 
proach to minimize the impact of these problems. We 
adopted severe criteria in the shape analysis with respect 
to the Paczyhski fit (stringent cut for both the and 
the Durbin- Watson test) and, furthermore, we eliminated 
candidates with both a long timescale > 40 days) 

and a red colour {{R ~ I)c > 1), since these most likely 
originate from variable stars. 

As we had already noted, this analysis suffered from 
the intrinsic limitation of an insufficient baseline, so that 
the variable-source issue could not be conclusively tested. 

Here we take advantage of the opportunity afforded by 
access to reduced INT data to extend our selected light 
curves. As we show, this will be a key element of our anal- 
ysis. 

3. 2. INT extension of MDM light curves 

The INT fields almost completely cover the MDM fields 
(except for a sizable fraction of the Control field far from 
the M31 bulge and a narrow band in the Target field). The 
good sampling during the two INT seasons 1999-2000 of 
data acquired in Sloan r' and i' bands^, together with 
their high quality, allows a straightforward and quantita- 
tive analysis. 

As a first step we have calculated the astrometric^ 
transformations between the different fields (data from 4 
different INT CCDs are needed to cover the two MDM 
fields). We found that a global transformation over one 
entire field led to systematic errors as large as 5". Hence, 
we use our large samples of 1194 and 845 references stars 
in the Target and Control field, respectively, to make a 
linear (i.e., 6-parameter) local transformation for each se- 
lected pixel. We demand a minimum of 24 reference stars, 
which are usually found within a square of ~ 2' — 3' about 
the selected pixel. We thereby reach an average astromet- 
ric precision of 0."2. To this we add the error in locat- 
ing the center of the variation, so that on average we get 
a ~ 0."5. 

As a second step we have determined from the cross- 
analysis of about 30 selected resolved stars in each field, a 
"colour equation" that allows us to align the MDM {R and 
/) and INT (r' and i') instrumental magnitudes closely 
enough to permit a quantitative comparison of the relative 
flux variations in each light curve. 

The uniqueness condition for a given variation requires 
that the extension of the light curve into the INT data 
should remain flat. However, the analysis of the PA-Nl 
event l|Auriere et al. 20T)T)l has shown that a microlensing 
light curve may be contaminated by a neighbouring vari- 
able star. Therefore, care must be taken before rejecting a 
microlensing event because of a second bump, and it is es- 
sential to test the following issues: the relative astrometric 

^ For this analysis we do not use the available g' band data. 
The INT pixel size is 0."33 versus 0."50 for MDM pixels. 



position of the two bumps and the similarity of the shape 
of the two deviations. 

Given a candidate microlensing light curve in the 
MDM data, we locate the corresponding INT light curve, 
and then calculate the estimator L in a square of 7 x 7 
pixels around the central pixel to check whether there is a 
variation, and where it is located. If we detect a variation 
in the INT data, we calculate its amplitude and width 
with a Paczyhski fit in order to compare quantitatively 
the two variations. 

As a criterion of rejection for a given MDM candidate 
microlensing event for which we find a variation on its ex- 
tension into the INT data, i.e., a likely variable star, we 
demand that the position of the INT variation be compat- 
ible with the corresponding MDM variation within 3 a and 
that the relative widths and deviations of flux lie within 
6 cr as evaluated from the two independent Paczyhski fits® . 

4. Candidate selection 

4-.1. Target and Control fields: the first selection 

We now present the results in both MDM fields. Target 
and Control, of an analysis complemented where possible 
by the stability test using the INT data. 

The threshold value for the estimators of the signifi- 
cance of the bump are fixed as in (Paper I| {Q > 100 and 
Li > 100, but now we proceed as explained in ^2.2.2t ). 
The selection criteria are the same as in [Paperl] the se- 
lected light curves must have enough points along the 
bump (at least 4 points on both sides of the maximum, 
and 3 inside the interval to ±ti/2/2); for the Paczyhski fit 
we require that x^/dof < 1.5 and that the Durbin- Watson 
test be satisfied at the 10% level; we require that either 
h/2 < 40 days or {R - I)c < 1. 

Due to the different approach we follow in bump de- 
tection, in the Target field we find 1 more light curve in 
addition to the 5 we already reported in (Paper I| In the 
Control field we find 4 light curves compatible with mi- 
crolensing. 

We summarize the main physical characteristics of 
these 10 events in Fig. Q which shows Rmax vs ii/2 and 

(i? - I)c vs ti/2- 

As can be noted, with respect to the selected light 
curves in [Paper I[ we tend to lack short timescale can- 
didates among the new events and, moreover, 3 of the 
5 candidates lie at the boundary we have fixed for the 
colour-timescale compatibility with a microlensing signal. 

* We allow here a larger margin because we are aware that 
probing such an effect can be difficult, especially when deal- 
ing with variations showing both a long time width and a red 
colour, these being possibly due to red variable stars that do 
not necessarily show a strictly periodic and regular behaviour. 

^ We note that, in principle, for a cluster of pixels with Q > 
100, we can have Li < 100. 
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The INT data allow us to check the stability of all 5 
Target light curves and 2 of the 4 Control light curves. 
In Table ^ we summarize the results of this comparative 
analysis, relative astrometry and relative shape analysis. 

As can be seen all 8 MDM events show compatible 
variations on their INT extensions. (In particular, we note 
that the positions of the two bumps are always compatible 
at 1 (T level). As a result, all the checked light curves are 
rejected as possible microlensing candidates. As an exam- 
ple, we show MDM light curves T4 and T5 (respectively, 
the shortest and the brightest flux variations detected) 
together with their INT extensions (Fig. 12)). 

This analysis shows how the sample of microlensing 
candidates derived from this selection is strongly contam- 
inated by variable stars. The two additional Control light 
curves for which there are no INT extensions have parame- 
ters: ti/2 = 51±5, 64±6 days and {R-I)c = 0.9±0.2, i.e., 
they are relatively long candidates. It then follows that, 
without INT vetting, it is no longer reasonable to look at 
them as viable microlensing candidates. Since, also look- 
ing at the results of the previous analysis, we suspect them 
to be red variable stars but cannot prove this without an 
extended baseline, we must exclude, when searching for 
such long variations, the non-INT portions of the MDM 
fields from the analysis. 

4-.2. Search for short duration events 

The analysis carried out in the previous section shows that 
variable sources can mimic the Paczyhski shape quite well, 
even with data available in two bands. It seems therefore 
appropriate in the search for viable microlensing candi- 
dates to relax the criteria introduced to characterize the 
shape of the variation and, on the other hand, to restrict 
the allowed space of physical parameters. We remark that 
the longer baseline now available make this approach vi- 
able. This search is also motivated by the lack of self- 
lensing events, which we would expect to find. 

We then proceed to a new selection with the same 
threshold values for L and Q. For the temporal sampling, 
we require at least 3 points inside the interval to i ti/2/2, 
with at least one on each side of ^o- We then impose the 
following conditions: 

- x^/dof < 5; 

- DW test at 2% significance level; 

- ti/2 < 20 days^*^. 

This search yields 8 additional light curves (4 in each 
field) which have durations ti/2 € (13,20) days and flux 
deviations Rmax G (21.0,22.8). The stability test on the 
INT data, as outlined in H3.2I allows us to reject 5 of these 
variations as microlensing candidates. We are left with 

We note that these short timescales are consistent with 
what we expect from the Monte Carlo simulations discussed in 
[Paper l| 



three light curves, all lying in the Control field. We label 
these C3, C4 and C5, the last having no INT extension. 
The issue of the stability on the INT extension of the two 
remaining candidates deserves some additional comments. 

MDM-C3: An INT variation is detected only in the i' 
light curve. Its position is compatible within 1 cr to that of 
the MDM variation. We note that the observed i' variation 
is significantly smaller than that observed on the MDM 
/ light curve and, as mentioned, the r' light curve is flat. 
Moreover, the time width of the MDM variation is actually 
quite short. This analysis indicates that, although sitting 
in the same position, the two variations may be due to 
different sources. At the same time we are aware that the 
INT data do not allow us to fully characterize the shape 
of the bump. We thus consider the stability test for this 
candidate to be inconclusive. 

MDM-C4: We detect a variation on the INT data that 
we localize within 1 a from the corresponding MDM varia- 
tion. We evaluate the variations of flux in r' and i' as being 
compatible within 3 and 1 a respectively with the corre- 
sponding variations on the MDM light curve. However, 
the observed time width of the INT variation, ti/2 ^ 70 
days, is significantly larger than the evaluated time width 
of the MDM variation, ti/2 ^ 13 days. As is the case for 
C3, the shape analysis may indicate that we are observing 
on the same light curve a variable star (a likely long pe- 
riod red variable) and a microlensing event. Our data are, 
however, insufficient to confirm or to reject this hypothe- 
sis. In this case too we then consider the stability test for 
this candidate to be inconclusive. 

In Table El we report the main characteristics of these 
light curves, position and physical parameters as evaluated 
from the Paczyhski fit. We note that, lacking any infor- 
mation from other sources (e.g. HST images) to measure 
the flux of the unamplified source, and since, on the ba- 
sis of the fit alone, the data do not allow us to break the 
parameter degeneracy, we have no way to get any reliable 
information on the physical parameter tg, so that also no 
reliable estimate of the mass of the lens is possible. 

As can be noted by looking at the light curves (Figs. 
13 El and EJ, the sampling does not allow us to test conclu- 
sively the symmetry of the bump. Regarding the achro- 
maticity, the ratio of R and / deviations from the baseline 
for the points belonging to the bump is about constant, as 
expected for microlensing events. We note that the MDM- 
C3 light curve, whose variation is quite short and very red, 
occupies a peculiar position in the (f? — I)c~ti/2 param- 
eter space (Fig.^ upper left of right panel). 

Even if the available data do not allow us to draw firm 
conclusions on the nature of the selected flux variations, 
their projected distance from the centre of M31 (see Table 
^ make them marginally consistent with the self-lensing 
events that we expect according to the Monte Carlo sim- 
ulations discussed in |Paper"l| 
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4.3. Comparison with POINT-AGAPE results 

In the framework of microlensing searches towards M31, 
the most promising results have been so far reported 
by the POINT- AGAPE collaboration. In particular, they 
have identified 4 high S/N (likely to be microlensing) 
light curves with high flux variation, Rmax < 21, and 
short timescale {ti/2 < 25 days, with 3 that are as 
small as ti/2 ~ 2 days), two of which lie quite near 
the M31 bulge, where, however, data still do not allow 
one to answer unambiguously the question of whether 
they are due to self-lensing effects or to MACHOs 
IjPauLin-Henriksson et al. 2003|l . We note that we lack, in 
the MDM data, such short duration and highly ampli- 
fied events. We conjecture that this apparent discrepancy 
can plausibly be explained by the relative quality of the 
data sets In particular, relative to the MDM data set, the 
INT data cover a much (> 3 times) larger portion of the 
sky around the M31 bulge, with more than the double 
the nights of observation. In particular, in the MDM Con- 
trol field, with at best one point every 3 nights, we can 
obviously not characterize events of such short duration. 
On the other hand, we point out that while the POINT- 
AGAPE analysis has been carried out with a threshold 
on the magnitude at maximum of the detected variations 
{Rmax < 21), in the analysis presented here we do not 
adopt such a threshold (we recall, see [Paper I[ that we 
are sensitive to variations down to around Rmax ~ 23). 

5. Discussion and perspectives 

We have reported on the status of the analysis of the 
MDM data carried out in the framework of a search for 
microlensing events towards M31. As a result of a full 
analysis of 2 years of data, in which we have also taken 
advantage of the longer baseline made available by the use 
of some INT data, we have first excluded as viable all the 
microlensing candidates previously reported in [Paper I| 
Second, we have reported the selection of three more can- 
didates compatible with a Paczyhski light curve, which 
may possibly be due to self lensing. On examining their 
extensions into the INT data, we find a second variation 
in both cases for which such data are available. However, 
these variations are strikingly different from and much 
longer than the short-duration events seen in the MDM 
data. Hence, it is plausible that each of these second vari- 
ations is due to an unrelated variable star superimposed 
on a true microlensing signal. The available data do not 
permit us to resolve this question. 

We have discussed the issue of variable-star contam- 
ination of the signal: we note that this is a significantly 
more serious problem than in the cases of the search for 
microlensing events towards the Magellanic clouds and the 
Galactic bulge in which one monitors the flux variations 
of resolved sources of known type. 

The analysis reported in this paper shows once again 
the crucial role played by a frequent sampling of the data 



and the total baseline length, both of which are essential 
to getting meaningful results when extracting microlens- 
ing signals from the background of variable sources. In 
this light, the forthcoming full analysis of the detection 
efficiency, together with results from the new sets of data 
acquired in the 2001 and 2002 seasons at both 1.3m and 
2.4m MDM telescopes with a new 8K-CCD array (as well 
as the prospective 2003 campaign at the TTl telescope in 
the south of Italy, IBozza et al. 1999|l . should eventually 
give us the opportunity to draw firmer conclusions on the 
issue of the MACHO fraction in the halo of M31 galaxy. 
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id 


(7(position) 




o-(A$b) 


cr(A<I'j) 


Tl 


1 


1 


2 


2 


T2 


1 


1 


1 


2 


T3 


1 


4 


4 


3 


T4 


1 


4 


5 


3 


T5 


1 


3 


2 


3 


T6 


1 


2 


1 


1 


CI 


1 


1 


2 


4 


C2 


1 


2 


2 


2 



Table 1. Results of the stability analysis on the INT data extension of the selected MDM microlensing candidates. 
For each MDM light curve (T and C stand respectively for Target, where we use the same numeration used in |Paper"I| 
and Control) for which wc find an astrometrically compatible INT bump, we report the number of standard deviation 
within which the position and the bump parameters (duration and flux deviation in both bands as calculated from 
independent Paczyhski fits on each of the MDM and INT light curves) , of the two variations are compatible. 
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d 




8'01" 


10'53" 


9'45" 


xVdof 




2.90 


2.32 


2.47 


dwR 




2.48 


2.01 


1.78 


dwi 




1.62 


2.23 


1.92 


ti/2 (days) 




16 ± 2 


13 ±2 


14 ±2 


Rmax 




21.0 ±0.1 


21.3 ± 0.1 


21.8 ± 0.1 


R-I 




2.2 ±0.2 


1.1 ±0.2 


0.5 ±0.2 


INT data 




inconclusive 


inconclusive 


no extension 



Table 2. Main characteristics of the three selected short-event light curves; d is the projected distance from the centre 
of M31. 
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Fig. 1. Physical characteristics for the selected light curves in i j4.1l (circles). Filled symbols represent candidates 
previously reported in Target field. Squares represent the result of the analysis presented in M.2I In the right panel, 
the lines at ti/2 = 40 days and {R — I)c = 1 indicate our excluded region of this parameter space. 
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Fig. 2. MDM T4 (left) and T5 light curves together with their extensions into the INT data. On the y axis, flux is 
in ADU/s; on the x axis, time is in days, with the origin in J-2449624.5 (both data sets). For the MDM light curves 
the dashed line represent the result of the Paczyiiski fit. For the INT light curves, shown together with the solid 
line representing the baseline is a dashed line representing the level of the maximum deviation of flux found on the 
corresponding MDM light curve. 
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Fig. 3. The MDM-C3 light curve (top) and its corresponding INT extension. Notation is the same as in Fig. |21 
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Fig. 4. The MDM-C4 light curve (top) and its corresponding INT extension. Notation is the same as in Fig. [2| 
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Fig. 5. The MDM-C5 hght curve, for which no extension on INT data is available. Notation is the same as in Fig. 



